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In experiments searching for a non-zero electric dipole moment of trapped particles, frequency 
shifts correlated with an applied electric field can be interpreted as a false signal. One such effect, 
referred to as the geometric phase effect, is known to occur in a magnetic field that is non-perfectly 
homogeneous. The increase in sensitivity of experiments demands improved theoretical description 
of this effect. In the case of fast particles, like atoms at room temperature and low pressure, the 
validity of established theories was limited to a cylindrical confinement cell in a uniform gradient 
with cylindrical symmetry. We develop a more general theory valid for an arbitrary shape of the 
magnetic field as well as for arbitrary geometry of the confinement cell. Our improved theory is 
especially relevant for experiments measuring the neutron electric dipole moment with an atomic 
co-magnetometer. In this context, we have reproduced and extended earlier numerical studies of 
the geometric phase effect induced by localized magnetic impurities. 



I. INTRODUCTION 

A non-zero electric dipole moment (EDM) for a spin 
1/2 particle would violate both the parity (P) and time 
reversal (T) symmetry, it would also violate the CP sym- 
metry according to the CPT connection. Since CP vi- 
olation beyond the standard model (SM) is needed to 
explain the generation of the baryon asymmetry in the 
early universe, electric dipole moments are important in- 
puts for cosmological models and theories beyond the 
SM [T] . In particular, the neutron EDM has been shown 
to strongly constrain many such theories So far, 

the best measurement of the neutron EDM [4] , extracted 
from the precession frequency of trapped ultracold neu- 
trons, is compatible with zero. 

When designing modern experiments to search for a 
non zero electric dipole moment of the neutron, the con- 
trol of the magnetic field fluctuations is of great signif- 
icance. A successful strategy is the use of an atomic 
comagnetometer, where polarized atoms fill the neutron 
confinement volume. In this volume, a combination of 
a weak magnetic field Bq and strong electric field E is 
applied, either in a parallel or an antiparallel configura- 
tion. The precession frequency of the comagnetometer 
measures the magnetic field via the first order relation 
uj = jBq, where 7 is the gyromagnetic ratio of the atom. 
The comagnetometer measurement is then used to cor- 
rect the neutron spin-precession frequency for the mag- 
netic field fluctuations as in the performed experiment 
[4] using a mercury comagnetometer [5]. In turn, the 
comagnetometer could induce false effects in the extrac- 
tion of the neutron EDM. The most dangerous effects are 



'Electronic address: guillaume.pignol@lpsc.in2p3.fr 
^Electronic address: roccia@csnsm.in2p3.fr 



frequency shifts Sco(E) of the comagnetometer correlated 
with the direction and value of the electric field. In this 
case, a false EDM for the comagnetometer is generated 
by the E-odd component of the frequency shift 
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that transmits into a false neutron EDM Ad n = 
7n/ 7 ^Faisci where j n is the neutron gyromagnetic ratio. 
In particular, it has been recognized both theoretically 
and experimentally that such effects arise from geomet- 
ric phase [6l [7] : when a particle moves with a velocity 
v with respect to an electric field E, it sees effectively a 
motional magnetic field Exv/c 2 . As a result, each spin- 
polarized particle is submitted to a fluctuating field (that 
we will refer to as w/7) resulting from a combination of 
the magnetic field inhomogeneities (B x , B y , B z ) and the 
motional field. Within this article we assume that the 
main magnetic field Bo — u>o /-y and the electric field are 
aligned with the z axis. The fluctuating transverse com- 
ponents Of IjJ 
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are known to induce a geometric phase shift in the spin- 
precession frequency (in the case of a rotating transverse 
field, the frequency shift is also know as the Ramsey- 
Bloch-Sicgcrt shift). 

The particular case of a uniform magnetic field gradi- 
ent with cylindrical symmetry 
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has been considered in [7]. In the low field regime, the 
resulting iS-odd frequency shift has been derived using 
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an elementary approach: 
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The calculation is done assuming a cylindrical trap of 
diameter D with the axis of the cylinder aligned with 
the z axis. The approach used in [7] is based on solving 
explicitly the Bloch equations for specific particle trajec- 
tories and could not be generalized to the case of arbi- 
trary shape of the transverse magnetic field. Then a more 
powerful theoretical treatment using the general theory 
of relaxation developed in [51 H] gave new insight to the 
problem, however still focusing on the case of a cylin- 
drical uniform field gradient. In this article we review 
and extend the general approach [5J |]5] , with the goal to 
apply it to the case of the mercury comagnetometer [5] 
used in the EDM experiment [10]. The magnetometer 
runs in the ballistic regime (gas collisions much less fre- 
quent than wall collisions) and in the nonadiabatic low 
field regime (the wall collision rate is much faster than 
the Larmor frequency) for which no theory was available 
beyond the case of a cylindrical cell in uniform magnetic 
field gradients with a cylindrical symetry. As a main re- 
sult of the present work, we establish the expression for 
the E-odd frequency shift valid for arbitrary shape of the 
magnetic field and for arbitrary shape of the trap: 
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where the brackets refer to the volume average over the 
trap. The case of an arbitrary field shape has been con- 
sidered in 11 for the case of a rectangular cell, with 
focus on the case where the particles are diffusing by gas 
collisions (high density limit). 

We first recall the general approach for the frequency 
shift and establish the terms proportional to E, E 2 , B 2 
respectively, valid for any shape of the magnetic inho- 
mogeneity and arbitrary shape of the trap. Then we 
discuss implications for neutron EDM experiments using 
an atomic comagnetometer in the nonadiabatic regime 
(such as the mercury comagnetometer). Whenever the 
geometry of the trap has to be specified for numerical 
studies, we will consider that of the experiment [10) (the 
OILL trap), a cylindrical trap whose axis is aligned with 
the main magnetic field, with a diameter D — 47 cm and 
a height H — 12 cm. 



II. GENERAL RESULTS FOR FREQUENCY 
SHIFTS 



The frequency shift induced by a fluctuating transverse 
field is given by the Lamoreaux-Golub expression [8]: 

1 f°° 
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L Jo 
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where the sign of the second term has been corrected (see 
appendix [A] for a derivation of this expression). Here the 
brackets refer to the ensemble average of the quantity 
over all particles in the trap. In our case where an electric 
field is present, the transverse field is given by ^ and 
the general formula ^ should be expanded in powers of 
the electric field: 



Suj = 5ujb 2 + Suje 2 + Suj be- 
The ^-independent frequency shift reads: 
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d T sm(uj r) (B x (0)B x (T) + B y (0)B y (r)). (8) 



It does not induce a direct false electric dipole moment 
since it remains unchanged when the electric field is re- 
versed. However, it has to be taken into account in the 
EDM analysis procedure. The quadratic frequency shift 
reads: 
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It can induce a false electric dipole moment if the mag- 
nitude of the E field changes when it is reversed [12j . 

Finally, the most important frequency shift is the E- 
odd term involving a linear dependence in the electric 
field: 
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(10) 

It generates a false electric dipole moment for the comag- 
netometer: 
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As we have just shown, the frequency shifts 
Slob 2 ! 5uje2 , Sujbe involve Fourier transform (evaluated 
at the Larmor frequency) of correlation functions in- 
volving field and velocity components. In the abiabatic 
regime the Fourier transforms can be expanded in per- 
turbation series using integration by part: 

f °°dT cos(wot) {B x (0)v x (t)) = (12) 
[cos(ujot){B x (0)x(t))]™ + ua f °° dr sin(oj T) (B x (0)x(t)) 

where the second term vanishes in the nonadiabatic limit 
{ujqt c -C 1). Using the previous expansion ( 
we arrive at the expression ([5]), which trans 
false EDM: 
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valid for arbitrary field inhomogeneities in the nonadia- 
batic regime. 

In the presence of a cylindrical uniform gradient ^ 
in a cylindrical trap, one has to calculate the volume 
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to show that the general result ( 13 1 reduces to 
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which coincides with eq. (37) of [7]. Let us point out 
that, although the full correlation function depends on 
the details of the random particle motion (such as the 
degree of specularity of reflection at wall collisions), this 
dependence disappears in the nonadiabatic limit, where 
only the volume average of field components is involved. 
The same conclusion was arrived at in [7J in the case of 
a cylindrical uniform gradient, using numerical simula- 
tions. 



III. IMPLICATIONS OF THE LINEAR 
FREQUENCY SHIFT 

Let us now work out the implications of the linear fre- 
quency shift of the comagnetometer when extracting the 
neutron electric dipole moment. 
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FIG. 1: Enhancement of the geometric phase effect calcu- 
lated for the OILL trap for a magnetic dipole on axis as a 
function of the distan ce f rom the bottom plate of the trap. 
The analytical result ( |29| l is compared to the heuristic and 
numerical results [13] . 



A. Large scale magnetic inhomogeneities 

First we deal with large scale magnetic inhomogeneities 
arising for, for example, non-perfect symmetry of the coil 
generating Bq or non-pcrfect shielding of the external 
ambient field. It is useful to parametrize the field com- 
ponents by polynomials of x, y, z up to a certain order. 
These polynomials however cannot be arbitrary since the 
Maxwell equations must hold. A convenient way to find 
a minimal polynomial expansion is to consider a scalar 
magnetic potential <fr for the field: 



B(x,y,z) = grad 4>(x,y,z). 



(16) 



In order for B to satisfy Maxwell's equations, it is neces- 
sary and sufficient for <f> to satisfy the Laplace equation 
Acf> = 0. Then, the general first order inhomogeneities 
are described by the following second order harmonic 
polynomials: 



(x, V ,z) = B z+ G ^x 2 + ^y 2 -'- (G x + G y )z 2 
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The inhomogeneous part has in total 5 degrees of free- 
dom. It corresponds to the transverse fields: 
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This general uniform gradient can be used to fit a mag- 
netic map in the experiment |10j and is much more gen- 
eral than the sole cylindrical uniform gradient ([3| consid- 
ered in earlier studies [7j. Then, from eq. ( |13[ ) we arrive 
at the expression for the false EDM of the mercury atoms: 
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Given that in this case G x + G y = —{^§f), the previous 
expression is in fact no different from that of a cylindrical 
uniform gradient. 

It suggests the following strategy of data analysis, 
which in fact was used in [4]: a correction of the mea- 
sured EDM is applied to subtract the geometric phase 
shift of the mercury atoms, according to: 
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the volume average magnetic field gradient {^§f) being 
evaluated independently. This strategy is reinforced by 
the new result presented here, which shows that the cor- 
rection is valid for the general uniform gradient (JTsJ) and 
not only for the case of cylindrical uniform gradient Q. 



B. Small scale magnetic inhomogeneities 



We have shown that the correction (20) is efficient 



at removing the false geometric phase effect induced by 
large scale magnetic inhomogeneities, but as already no- 
ticed in [U [T3] , it could become inefficient in the presence 
of a localized magnetic inhomogeneity where the result- 
ing field gradients are not uniform over the trap volume. 
The correction (20) is justified by the fact that for a 



general uniform gradient ( 18 ) the following expression 
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holds in a cylindrical trap of radius R. In the general 
case, to account for localized magnetic field, we define 
the enhancement factor 1 + E as 
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This enhancement factor strongly depends on the prox- 
imity of the small scale inhomogeneity. Let us consider 
the case of a field induced by a dipole of magnetic mo- 
ment m = me z situated below the bottom electrode, on 
the axis of the cylindrical trap. The field reads, in cylin- 
drical coordinates centered at the dipole position: 
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where p = pqiti/Att is the strength of the dipole, d — 
\J ' z 2 + p 2 is the distance from the dipole. One can cal- 
culate the gradient 

dB z f9z 15z 3 \ . , 

ar^U-^J (25) 

and its volume average 
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where z d is the distance of the dipole from the bottom 
electrode. Now the integral of the left hand side of eq. 
([22} is 

z d +H 
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Then one can calculate analytically the enhancement fac- 
tor: 
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As shown in fig. [T] the analytical result derived here 
reproduces very well earlier numerical studies [13] . 

The improved theory can now predict the enhancement 
factor even when the dipole is at the surface of the cell, 
which corresponds to the worst case. We have calculated 
the systematic effect on the neutron EDM due to the 
presence of a dipole at the surface that remains after the 
correction (20). When the dipole is not situated at the 



center of the electrode, we did not find analytical expres- 



sions for the volume average (13), the volume integrals 
have been calculated with numerical methods. The re- 
sult is shown in fig. [2j as a function of the position of 
the dipole on the bottom electrode and for two orienta- 
tions of the dipole. We found that a dipole situated at 
the circumference is unfortunately particularly effective 
at generating a false EDM. This result can be used to set 
the requirements on the control of magnetic impurities 
in an experiment searching for the neutron EDM. 
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FIG. 2: False neutron EDM signal due to a dipole of strength 
p = 0.5 nTcm at the lower surface of the cell for different 
radial positions and different orientations. The left hand side 
plot shows the main contribution around the circumference of 
the cell. The right hand side plot is a zoom for smaller radii. 



IV. CONCLUSION 

In the context of the search for the neutron EDM with 
an atomic comagnetometer, most of the systematic ef- 
fects are related with the geometric phase shift. In the 
nonadiabatic limit which describes a comagnetometer at 
low pressure, an improved theory valid for any geometry 
of the trap and arbitrary shape of the magnetic field was 
developed. It permits to calculate the false EDM gener- 
ated by the geometric phase shift of mercury using the 
general formula ( 13 ) thus avoiding recoursing to heavy 



spin-tracking Monte-Carlo codes. 
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Appendix A: Spin relaxation and frequency shifts: 
an elementary derivation 

When an ensemble of spin-polarized particles evolve 
in a weakly fluctuating magnetic field, two important 
phenomena occur: spin-relaxation and shift of the res- 
onance line. The situation is particularly relevant for a 
gas of polarized particles evolving in a static, but inho- 
mogeneous magnetic field. In this case each particle sees 
effectively a time dependent magnetic field, which differ 
from particle to particle. Here, an elementary derivation 
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of the basic formulas is given for the spin 1/2 case. The 
classic derivations are usually done within the density 
matrix formalism, as a special case of the general theory 
of relaxation developed by Redfield and others (see for 
example [14] for such a derivation). Here we derive the 
same results using the classical equations of motion of 
single magnetic moments M with gyromagnetic ratio 7 
evolving in a weakly fluctuating magnetic field B: 
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We assume that the magnetic field can be written in the 
following form: 
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B(t) = B a e z + -uj{t) with (u>(t)) = 0. 
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It is indeed always possible to define the z axis to coincide 
with the direction of the average field. Here (X) denotes 
the ensemble average of the quantity X. We proceed 
and define the magnetic moments in the frame rotating 
around the z axis at the frequency ujq — 7 -Bo: 



m = R(t)M, 
with R(t) being the rotation matrix: 
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Then the magnetic moment of a given particle seen in 
the rotating frame evolves as: 
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where b(t) = R(t)u)(t)/~{ is the fluctuating magnetic 
field in the rotating frame. We are interested in deriving 
a simple evolution equation of the average magnetization 
(m). We start with the ensemble average of eq. ( A5) 
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which can be successively integrated: 
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The first term vanishes because we assume the absence 
of correlation between the initial direction of m and the 
future magnetic field b for individual particles. The vari- 
ations of (m) will be considered as a perturbation, which 
is valid in the limit of small fluctuating fields (b 2 ) <§C By. 
The perturbative approach consists in substituting in the 
second term m(0) in place of m(t'). Then we arrive at 



the first order equation, by introducing a new variable 
T = t-t': 

d ( m > = 7 2 I d T (b(t) x (b(t - t) x m)) (A8) 
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= -r<m) 

where we define the relaxation matrix T. Using the triple 
product expansion, we obtain an expression of the re- 
laxation matrix involving autocorrelations of the compo- 
nents of the fluctuating field: 
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This matrix should be decomposed in the following form, 
with transparent meaning: 
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The second matrix in the right hand side does not lead 
to first order effects because it corresponds to oscillating 
terms, at the frequency 2wo, when transforming the equa- 
tion (A8) back to the laboratory frame. All important 



(A4) phenomenological consequences of the fluctuating field 
are encoded in the first matrix, where ri = ^- is the lon- 



gitudinal spin-relaxation rate, r 2 = tjt is the transverse 
spin-relaxation rate, and Sui represents the shift of the 
resonant frequency. 



Now, one can derive from (A9) and (A10) the relevant 



terms of the relaxation matrix, in terms of the autocorre- 
lation of the fluctuating field expressed in the laboratory 
frame. 
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To get the final expressions (All ) ( A12 1 (A13), further 



reasonable assumptions have been made. It has been 
assumed that the system is stationary in the statistical 
sense, that is, the correlations depends only on the time 
difference: 



(b i (t)b j (t-r)) = (b i (r)b j (0))- 
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Also, it has been assumed that the integrals in eq. ( A9 ) 



can be extended to infinite time, which is justified be- 
cause the correlations decay with a time much shorter 
than the observation time. When disregarding the tran- 
sient regime immediately after the time zero for an ob- 



servation, the relaxation matrix is time-independent. As 
a final remark, a similar calculation has been done in the 
density matrix formalism to obtain the frequency shift in 
jS], which agrees with (A13) after correcting [5] for a sign 
error. 
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